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Abstract 

We study synthetic aperture radar (SAR) imaging and motion estimation of complex scenes 
consisting of stationary and moving targets. We use the classic SAR setup with a single antenna 
emitting signals and receiving the echoes from the scene. The known motion estimation methods 
for such setups work only in simple cases, with one or a few targets in the same motion. We 
propose to extend the applicability of these methods to complex scenes, by complementing them 
with a data pre-processing step intended to separate the echoes from the stationary targets and 
the moving ones. We present two approaches. The first is an iteration designed to subtract the 
echoes from the stationary targets one by one. It estimates the location of each stationary target 
from a preliminary image, and then uses it to define a filter that removes its echo from the data. 
The second approach is based on the robust principle component analysis (PCA) method. The 
key observation is that with appropriate pre-processing and windowing, the discrete samples of 
the stationary target echoes form a low rank matrix, whereas the samples of a few moving target 
echoes form a high rank sparse matrix. The robust PCA method is designed to separate the 
low rank from the sparse part, and thus can be used for the SAR data separation. We present 
a brief analysis of the two methods and explain how they can be combined to improve the data 
separation for extended and complex imaging scenes. We also assess the performance of the 
methods with extensive numerical simulations. 



1 Introduction 

Synthetic aperture radar (SAR) is an important technology capable of computing high resolution 

x 

images of a scene J on the ground, using measurements made from an antenna mounted on a 
platform flying over it (HI HI E)J. The antenna periodically emits a signal f(t) and records the 
echoes, the data that are processed to form an image. The problem setup is illustrated in Figure [T] 
The data D(s,t) are indexed by the slow-time s and the fast-time t. The slow-time parameterizes 
the antenna position r(s) at the moment it emits the signal. The fast-time parameterizes the time 
between consecutive illuminations at s and s + As, so that < t < As. 
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Figure 1: Setup for synthetic aperture imaging. 

The recordings D{s 1 t) are approximately and up to a multiplicative factor a linear superposition 
of the probing signals f(t) time-delayed by the round-trip travel time between the antenna position 
r(s) and the locations p of the scatterers (targets) in J . Assuming that the waves propagate at 
the constant speed of light c, the round-trip travel time is given by 

r(s,p) = 2\r(s)-p\/c. (1.1) 

To form an image, the data are first convolved with the time reversed emitted signal /(£), a 
process known as match-filtering or pulse compression. Then, they are back-propagated to search 
points p x E J X using the travel times r(s,p x ). The image is obtained by superposing the pulse 
compressed and back-propagated data over a slow time interval defining the synthetic aperture a, 
the distance traveled by the antenna. 

SAR systems that cover long apertures and emit broad-band high frequency signals give very 
good resolution of images of stationary scenes. The resolution can be of the order of ten centimeters 
at ranges of ten kilometers away. However, the images are degraded when there is motion in the 
scene. Depending on how fast the targets move, they may appear blurred and displaced, or they 
may not be visible at all. The motion needs to be estimated, and then compensated in the image 
formation, in order to bring these targets in focus. This is a complicated task for complex imaging 
scenes, consisting of many stationary targets and moving ones. 

There are basically two approaches to motion estimation. The first determines the motion from 
the phase modulations of the return echoes [3 [TTJ [TUl HI EHl [26l [15l [TTl [20l [22l [12] , assuming either 
a single target in the scene, or many targets moving the same way. It does not work well in complex 
scenes. The estimation is sensitive to the presence of strong stationary targets and to targets that 
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have different motion. The second approach uses multiple receiver and/or transmitter antennas 
[T6t 125} [23] . It forms a collection of images with the echoes measured by each receiver-transmitter 
pair. Then, it extracts the target velocities from the phase variation of the images with respect to 
the receiver/transmitter offsets. 

We are interested in the first approach that uses the classic SAR setup with a single antenna. 
To extend its applicability to complex scenes, we complement it with a data pre-processing step 
intended to separate the echoes from the stationary targets and the moving ones. A successful 
separation allows motion estimation to be carried with the echoes from the moving targets alone, 
using the algorithms in mMDHHEMHliP- 

We present two data separation methods. The first extends the ideas in |3J |5], developed in 
a different context of imaging in scattering layered media. The method seeks to subtract the 
stationary target echoes one by one, using a combination of travel time transformations and a 
linear annihilation filter. Each travel time transformation is relative to one stationary target at a 
time, whose location can be estimated from a preliminary image. After the transformation, the 
echo from this target becomes essentially independent of the slow time, and so it can be annihilated 
with a linear filter such as a difference operator in s. The travel time transformation is undone 
after the filtering, and the method moves on to the next stationary target. 

The second method is based on the robust principle component analysis (PCA) method [7J. We 
have shown in [2] with analysis and numerical simulations that after appropriate pre-processing and 
windowing, the samples of the echoes from stationary targets form a low rank matrix. Moreover, 
the samples of the echoes from a few moving targets form a high rank but sparse matrix. Thus, the 
echoes can be separated with the robust PCA method, which is designed to decompose a matrix 
into its low rank and sparse parts [7J. We review briefly the results in [2]. We also explain that 
in practice it may be useful to use a combination of the two approaches in order to achieve better 
data separation for extended and complex imaging scenes. 

The paper is organized as follows. We begin in section [2] with a brief description of basic SAR 
data processing and image formation. The data separation methods are described in section |3j We 
analyze them in section [4] for simple imaging scenes. Then we show with numerical simulations in 
section [5] that the methods can handle complex scenes. We end with a summary in section [6| 

2 Preliminaries 

We begin with two data pre-processing steps: pulse compression and range compression. They are 
commonly used in SAR, but in addition, they are an important first step in our data separation 
methods. We also describe the generic SAR image formation. The goal of the paper is to explain 
how to complement the image formation with data separation and filtering in order to estimate 
motion and focus images of complex scenes. 
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2.1 Pulse and range compression 

Typically, the antenna emits signals that consist of a base-band waveform modulated by a 

carrier frequency vq = uj / (2tv) 

f(t) = cos(u t)f B (t). (2.1) 

The Fourier transform of the signal is 

— f°° 1 — i 

f(u>)= dtf(t)e Mt = - f B (e> + u> ) + f B (u>-u> ) , (2.2) 

J-OO Z 1 J 

where the support of is the interval [—tvB^tvB], and B is the bandwidth. Ideally, the 

waveform /b(£) should be a pulse of short duration, so that travel times can be easily estimated. 
It should also carry enough power so that the received echoes are stronger than the ambient noise. 
However, antennas have instantaneous power limitations and they can transmit large net power 
only by spreading it over long signals, like chirps. Consequently, the measured echoes are long and 
it is impossible to read travel times directly from them. This difficulty is overcome by compressing 
the echoes, as if they were due to short emitted pulses. 

Pulse compression is achieved by convolving (match filtering) the data with the complex con- 
jugate of the time reversed emitted signal 

D p (s,t) = J dt'D{s,t')f{t f -t). (2.3) 

Mathematically, this is equivalent to the antenna emitting the signal 

f p (t) = J dt'f(t')W~t), (2.4) 

and receiving the echoes D p (s, t). The signal f p (t) is a pulse of duration of the order l/£>, as shown 
for example in [18J. 

Range compression is another important data pre-processing step. It consists of shifting the 
fast time t in the argument of D p {s, t) by the travel time to a reference point p Q E J , 

D r (s,t') = D p (s,t' + r(s,p )). (2.5) 

Here t' is the shifted fast-time, satisfying t' + r(s, p Q ) E [0, As]. 

Range compression is beneficial for numerical computations and it is essential for our data 
separation methods because it removes partially the dependence of D p (s,t) with respect to the 
slow time. We illustrate this in Figure [2j where we compare the pulse compressed echo from a 
stationary target before and after the range compression. The target is at location p = (0,5,0)m 
and the reference point is at p Q = (0,0, 0)m. The SAR platform moves at speed 70m/s along a 
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Figure 2: Left: Pulse compressed data D p (s,t) = f(t — r(s,p)) from a target at p = (0, 5,0)m, 
recorded by a SAR antenna moving at speed V = 70 m/s along a linear aperture of 400 m, at range 
L = 10 km from p. Right: Range-compressed data D r {s 1 t r ) with respect to p Q = (0,0, 0)m. 

linear aperture a = 400m at range L = 10 km from p. The imaging scene lies on a flat surface at 
zero altitude. We show in the left plot in Figure [2] the amplitude of the pulse compressed echo as 
a function of s and t. Note that its peak lies on the curve denned by 




This curve is a hyperbola because the aperture is linear and the platform moves at constant speed. 
The right plot in Figure [2] shows the amplitude of the range and pulse compressed echo. Its peak 
lies on the curve 

ct! 

— = \r(s)-p\ - |r(» -pj, 

which looks almost as a vertical line segment, because the dependence on the slow time has been 
approximately removed by the range compression. 

We work henceforth with the range and pulse compressed data and drop the prime from the 
shifted fast time t' to simplify the notation. Borrowing terminology from geophysics, we call D r (s, t) 
the data traces or simply the traces. 

2.2 Image formation 

A SAR image is formed by superposing the traces over the aperture and back propagating them to 
the imaging points p using travel times 

n/2 

1(P X )= J2 D r { S 3^{ S ^P X ) -<S h p )). (2.6) 

j=-n/2 
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The slow time samples sj discretize the interval [— S(a), S(a)] which defines the aperture a along the 
trajectory of the SAR platform. The sampling As is uniform and equals the time interval between 
tho consecutive signal emissions from the antenna. To simplify the notation, we assume that the 
origin of the slow time corresponds to the center of the aperture. Thus, n is an even integer and 

2S{a) = nAs. 



Equation (2.6) is a simplification of the imaging function used in practice, which includes a 
weight factor that accounts for geometrical spreading. The weight makes a difference only when the 
aperture is large. We neglect it in this paper because we focus attention on motion estimation, which 
is necessarily done over small successive sub-apertures. Targets may have complicated trajectories 
over long data acquisition times, but we can approximate their motion by a uniform translation 
over short durations 2S(a) defining small sub-apertures a. The idea is that the speed of translation 
can be estimated sub-aperture by sub-aperture, and then it can be compensated in the image 
formation. By compensation we mean that if we seek to image a moving target and we estimate 
its speed to be u, we can form an image as 

n/2 

X aO X )= Yl D r ( Sj ,r ( Sj ,p X + Sj uj -r( Sj ,p )y (2.7) 

j=-n/2 

We repeat the process for successive sub-apertures, and assemble the final image by superposing the 
results. This superposition requires proper weighting to compensate geometrical spreading effects 
over long SAR platform trajectories. 

Motion estimation and its compensation in imaging works for simple scenes with one or a few 
targets that move at approximately the same speed. To our knowledge, there is no approach that 
estimates motion in complex scenes with a single antenna. Moreover, even if a target speed were 



estimated, the compensation in (2.7) would bring that target in focus but it would blur the others. 
To address the difficulty of imaging with motion estimation in complex scenes, we propose a data 
pre-processing step for separating the traces from the stationary targets and the moving ones. If 
we could achieve such separation, we could carry the motion estimation on the data traces from 
the moving targets alone. Moreover, we could image the stationary and moving targets separately, 
and then combine the results to obtain an image of the complex scene. 



3 Data separation 

We present two complementary approaches for separating traces from stationary and moving tar- 
gets. They are analyzed in section |4j in simple setups. Here we describe how they work and 
illustrate the ideas with numerical simulations, for a complex scene named hereafter "Scene 1". 



6 





Figure 3: Left: The image of the complex Scene 1, with twenty stationary targets and two moving 
targets. Right: Image of the same stationary scene, and no moving targets. The effect of the 
moving targets in the left image amounts to two faint streaks. 

It consists of twenty stationary targets and two moving targets with speeds of 28m/s and 14m/s, 
respectively. We refer to section [5] for details on the setup of the numerical simulations. The traces 
are displayed in the left plot of Figure |6j 

3.1 Travel time transformations and data separation 

The basic idea of the first data separation approach is that if we have an estimate p e of the location 
p of a stationary target, we can shift the time t by the travel time r(s, p e ) to remove approximately 
the dependence on the slow time s of the trace from that target. Then, we can subtract the trace 
from the data by exploiting its independence of s. To determine p e , we need a preliminary image 
of the scene. The assumption is that the stationary targets dominate the scene, and therefore they 
are visible in the preliminary image even in the presence of moving targets. See Figure [3] for an 
illustration. 

-+e 

Let T:j_ be the map taking the range compressed data D r (s,t) to the data in the transformed 
coordinates (s, t + Ar(s, p 6 )) 

[Tf D r ](s, t) = D r (s, t + Ar(s, p e )), (3.1) 

where 

AT(s,p e ) = T(s,p e )-T(s,p ). (3.2) 

Recall that p Q is the reference point in the range compression. To illustrate the effect of the map 
(3.1), we show on the left in Figure [4] the trace from a target at location p ^ p Q . The plot on the 



right shows the trace after the travel time transformation (3.1), with the ideal estimate p e = p 



The transformation makes the trace independent of s, and thus it can be annihilated by a difference 
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Figure 4: Left: The trace \D r (s,t)\ from one stationary target. Right: the traces after the travel 
time transformation relative to the exact location of the target. Note that there is no dependence 
on s after the transformation. 



operator in s, 



D s [Tf D r ](s,t) = {[Tf D r ](s + As,*) - [Tf D r ](s,t)} « J^[t£ D r ](s,t). (3.3) 

This is an extension of the ideas proposed in [H |5] for annihilating echoes from a layered scattering 
medium. Other annihilation operators than ( |3.3[ ) may be used, as well. In practice, where noise 
plays a role, they should be complemented with a smoothing process. We refer to [U IS] for a 
detailed analysis of annihilator filters. For the simulations in this paper the annihilator D 5 works 
well, and there is no noise added to the data. 

After taking the derivative, we map the traces back to the (s, t) coordinates using the travel 
time transformation T_ . The data filter is a linear operator denoted by Q p , given by composition 



of T£ , D s and 



0? D r 



(M) 



oT?. D r 



(M) = [B s D r (s,t f + AT(s,p e ))] t , =t _ Ar{s , ey 



(3.4) 



We demonstrate in Figure [5] the travel-time transformations and annihilation steps for a scene with 
one stationary and one moving target. Note that the trace from the stationary target located at p 
is completely removed by the mapping , because p e = p. 

If there are N stationary targets in the scene, we can subtract their traces from the data one 
by one, by iterating the procedure above. Let p^ be the estimated locations of the targets, with 



j = 1, . . . , N. The data filter Q is the composition of the linear operators (3.4) 



[QD r ] (s,t) 



QPn o QPn-i o . . . o QP* D r 



(M)- 



(3.5) 



The estimates p e , are not exact in practice, but the filter (3.4) is robust with respect to the 



estimation errors. We show this in section |1J where we present a brief analysis of the filters. In 
fact, our simulations show that it is typical that one application of QP removes at once the traces 
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Figure 5: Left: The traces D r (s, t) for a scene with one stationary target at location p and a moving 
one at 28m/s. Middle: The traces after the travel time transformation [T^_D r ](s,t). Right: The 
filtered traces [Q? e D r ]{s,t). 



from all the stationary targets located near p e . Thus, the needed number of iterations in (3.5) may 
be much smaller than the number of stationary targets in practice. 

Because antennas measure discrete samples of the data D r (s,t), we need to use interpolation 
when making the travel-time transformations. We do so by implementing the travel-time transfor- 
mations as phase modulations in the frequency domain 



where 



( s ,t) = 1-Jdu D r (s,oj)e-^ Ar ^\ (3.6) 
D r (s,tj) = J dt e iujt D r (s,t). 



We implement (3.6) in the discrete setting using the fast Fourier transform. 



3.2 Data separation with robust PCA 

We introduced recently in [2] a method for separating the traces from stationary and moving targets 
using robust principal component analysis (PCA). We review here the data separation method, and 
give some details of its analysis in section [I| 

The robust PCA method was introduced and analyzed in [7]. It decomposes a matrix M into 
a low rank matrix C and a sparse matrix S. In our context, the entries in M are the samples of 
the traces 

M j£ = D r (s j ,t £ ), (3.7) 

where 

sj = jAs, j = -n/2, ...,n/2, 2S(a) = nAs, (3.8) 

and 

U = £As, t = -ra/2, ...,m/2, As = mAt. (3.9) 
The sampling is at uniform intervals As and At, determined in practice by the SAR antenna 
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Figure 6: Left: The matrix M of the sampled traces from Scene 1. Middle: The matrix C consisting 
of the traces from the stationary targets in Scene 1. Right: The matrix S consisting of the traces 
from the moving targets in Scene 1. 



system. To simplify the notation, we take the origin of s at the center of the aperture, and the 
origin of t in the middle of the interval As between two consecutive signal emissions. 

To explain why we can use robust PCA for the data separation, consider the matrix M. of 
sampled traces for Scene 1. We display it in the left plot of Figure |6j The traces that appear 
almost vertical correspond to the stationary targets. This is because these targets are not too far 
from the reference point p G , and the range compression removes most of their s dependence. We 
expect that these traces form a low rank matrix C. We display them in the middle plot of Figure [6j 
The range compression does not remove the slow time dependence of the traces from the moving 
targets, so they appear sloped in Figure [6| They form a high rank but sparse matrix <S, shown in 
the right plot of Figure |6j 

The robust PCA method is designed to separate C from S in M = C + S. It does so by solving 
an optimization problem called principle component pursuit [7] 

min l|£|L + 77||S||i, subject to £ + S = M. (3.10) 

£,<SGM n + 1Xm + 1 

Here ||£||* is the nuclear norm of £, the sum of the singular values of £, ||5||i is the matrix 1-norm 
of S and 

rj= , 1 (3.11) 

y max{n + 1, m + 1} 



The optimization (3.10) can be applied to any matrix A4. But it is shown in [7] that if M is a 
matrix with a low rank part C and a high rank sparse part 5, then the algorithm recovers exactly 
jC and S. We refer to [7j for sufficient (not necessary) conditions satisfied by so that it can be 
succesfully decomposed by principle component pursuit. 

Robust PCA is well suited for our purpose, but it cannot be applied as a black box to the matrix 
of traces. We need to complement it with properly calibrated data windowing. To illustrate why, 
we present two examples. 
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Figure 7: Separation of traces from the scene in Example 1. From left to right we show the matrix 
A4 of traces, the low rank £, and sparse S components obtained with robust PCA. Top row: robust 
PC A applied to the entire data matrix M. Bottom row: robust PCA applied to a windowed M. In 
each plot the matrices are normalized by the largest value of \M\. The sparse components shown 
in the right column are plotted in (decibel) dB scale to make the contrast more visible. 

Example 1: The first example considers a scene with one moving target and seven stationary ones. 
The moving target is at location (0, 0, 0)m at time 5 = 0, corresponding to the center of the aperture, 
and velocity u = 28/v^2(l, 1, 0)m/s. The stationary scatterers are located at (0, 0, 0)m, (±5, 0, 0)m, 
(0, ±5, 0)m, and (±10, 0, 0)m. in the imaging region. The reference point is p Q — (0, 0, 0)m. 

When we apply the robust PCA method to the matrix of traces shown in the top left plot in 
Figure [7J we obtain the "low rank" and "sparse" parts shown in the top middle and top right plots. 
We note from the right plot that the "sparse" component contains the trace from the moving target 
but also remnants of the traces from the stationary targets. The data separation is not successful 
because the matrix Ai is sparse to begin with. 

However, if we apply robust PCA to the matrix of traces in a smaller fast-time window, we get 
a successful separation. This is shown in the bottom row of Figure [7| The windowed A4 has a 
much larger number of non-zero entries relative to its size, and robust PCA separates the traces 
from the moving target and the stationary ones. 

Example 2: The second example considers a scene with thirty stationary scatterers and the same 
moving target as in example 1. The matrix M of traces is displayed in Figure |5J The results with 
robust PCA applied to the entire M are in the top row of Figure |9| Again, we see that the "sparse" 
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Figure 8: Matrix M. of traces from the scene in Example 2. 
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Figure 9: Robust PC A separation of traces from a scene with thirty stationary targets and one 
moving target. The top row of plots shows the results with robust PC A applied to the entire 
data matrix in Figure |8j The low rank and sparse components are shown on the left and right, 
respectively. The bottom row of plots shows the results with robust PCA applied to successive 
small time windows of the data. The sparse components shown in the right column are plotted in 
dB scale to make the contrast more visible. 

part contains remnants of the traces from the stationary targets. This is not because M is sparse, 
as was the case in example 1, but because the traces from the stationary targets do not form a 
matrix of sufficiently low rank. When we window .M, we work with a few nearby traces at a time. 
These traces are similar to each other so they form a low rank matrix that can be separated with 
robust PCA. The final result is given by the concatenation of the matrices in each window. It is 
shown in the bottom row of Figure [9j 

These two examples show that windowing is a key step to successful data separation with robust 
PCA. To choose appropriate window sizes, it is necessary to understand how the rank and sparsity 
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of the data matrix M depend on the location and density of the stationary targets, and on the 
velocity of moving targets. This is addressed in the analysis of [2], which we review in section 4.4 



3.3 Discussion 

We use numerical simulations for Scene 1 to compare the two methods described in sections 3.1 



and 3.2 The results are in Figures 10 and 11 
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Figure 10: Annihilation with the filter (3.5) of the stationary target traces from Scene 1. Top: the 



matrix of traces. Bottom: the filtered traces. In each plot we show absolute values normalized by 
the largest value of \M\. The right column is the left column plotted in dB scale. 



The first method uses the filter Q defined in (3.5) to annihilate the traces from the stationary 



targets. We see from Figure [10| that this is accomplished, but the filter removes some parts of the 
traces from the moving targets as well. All our numerical simulations show that Q separates nicely 
the traces from the moving targets if the scene does not have too many stationary targets at almost 
the same range. Otherwise, it also removes parts of the traces from the moving targets. 

This 



The data separation with the robust PC A approach is better, as seen in Figure 11 



method is more sophisticated than the other, because it involves careful windowing of the traces, 
but it works better for complex scenes that are not too large. Robust PCA relies heavily on range 
compression removing most of the slow-time dependence of the traces from the stationary targets. 
But for extended scenes there is no single reference point p Q that accomplishes this task, and data 



separation with robust PCA deteriorates. The filter (3.5) still works for extended scenes, because 



the travel time transformations in (3.4) are relative to the location of the targets that we wish to 
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Figure 11: Trace separation with robust PC A for Scene 1. The low rank part C is shown on the 
left and the sparse part S is shown on the right. In each plot we show absolute values normalized 
by the largest value of \M\. The sparse part is plotted in dB scale. 

annihilate, and not some arbitrary reference point. 

In principle, we can use a combination of the two approaches to improve the data separation 
with robust PCA for extended scenes. The idea is to begin with a preliminary image and identify 
roughly the stationary targets to be separated from the scene. Then we can window the image 
around groups of such targets and work with one subscene at a time. It is important to observe here 
that if the aperture is not too large, as is the case in motion estimation, the image is approximately 
the Fourier transform of the data with a phase correction. See for example [3] for a detailed 
explanation of this fact. Thus, we can obtain approximately the traces from the smaller subscene 
by Fourier transforming the windowed image. To separate these traces, we can apply the travel 



time transformation (3.1) with a reference point p e in the subscene, followed by the robust PCA 



approach as described above. 

3.4 Velocity estimation and separation of traces from multiple moving targets 

Let us assume in this section that the traces from the stationary targets have been removed, and 
denote by the same D r (s,t) the traces from the moving targets. Since there may be multiple 
moving targets, we would like to separate these traces further, based on differences between the 
target velocities. We can do so by extending the annihilation filters described in section 3J_, as we 
explain next. 

We need a key observation, shown with analysis in section [4| The slope of the trace from a 
moving target is determined approximately by its speed along the unit vector 

a. = (3.12) 

l r (0) - Pol 

pointing from the SAR platform location r(0) at the middle of the aperture to the reference point 
p Q in the imaging scene. Moreover, the curvature of the trace is determined by the speed projected 
on the plane orthogonal to m . 
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Because we assume a flat imaging surface, the target speed u is a vector in the horizontal plane 



u (u, 0), u E 



(3.13) 



We parametrize it by the range velocity u and the cross-range velocity u±. The range velocity is 



u = u • m u • m G , 



(3.14) 



where m is the projection of m G on the imaging plane. The cross-range velocity is defined by 



u± = t ■ P G u = u • t — um • t(0), 



where P is the orthogonal projection 



(3.15) 



(3.16) 



and t is the unit vector tangent to the SAR platform trajectory, at the center r(0) of the aperture. 
Note that in many setups m and t are almost orthogonal, so u± is approximately the target speed 
along the vector t, the projection of t on the imaging plane. 



Next, we define the transformation mapping T" ' p , the generalization of (3.1) 

(s, t) = D r (s, t + Ar(s, p e + su e )) , 



(3.17) 



with At given in (3.2). Here p e is the estimated location of the moving target at the origin of the 



slow time, and u e is the estimated velocity. If the line segment 



p e + su e , 8€[-S(a),S(a)], 



is close to the target trajectory, then the transformation (3.17) removes approximately the depen- 



dence of its trace on the slow time. Thus, we can subtract it from the other traces using the same 



difference operator D s . The annihilation filter is the generalization of (|3.4|) 



(s,t) = 



T 



D r 



(s,t) 



s D r (s, t' + At(s, p e + su e ))] t , =f _ A 



T(s,p e +SU e ) ' 



(3.18) 



where the mapping T_ ,p undoes the travel time transformation. Additionally we can apply 



robust PCA to 



Tj 



D r 



This should result in the low rank matrix containing the trace 
corresponding to the moving target with velocity p e and the sparse matrix containing the rest of 
the moving target traces. 



The implementation of (3.17) and (3.18) requires the estimates p e and u e . There are two 
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scenarios. The first assumes that we know p e from tracking the target in previous sub-apertures, 
and seeks only u e . The second seeks both p e and u e . In either case, the estimation begins with 
that of the range velocity, which controls the slope of the trace. Then, we can estimate p e if we 
need to, and the cross-range speed u\. 

We estimate u e as the maximizer of the objective function 



max 

-m/2<£<m/2 

3 



n/2 

E \[t 

=-n/2 



D r 



(sj,ti) , u = u(m o ,0). 



(3.19) 



The idea is that when u e is close to the range speed of the target, the transformation T"'^ makes 
the trace an approximate line segment, and when we sum it over s, we get a peak around some 



fast time for \l\ < m/2. If we do not know p e , we can set it equal to p Q in (3.19). Then, we can 
form an image Z^e(p x ) using (2.7) with u e = tx e (m o ,0), and estimate p e from there. The motion 



compensation in ( |2.7[ ) brings the moving target approximately in focus, and this is why we expect 
to be able to estimate its location from the image. 

Once we have estimated u e and p e , we seek the cross-range speed u\ as the minimizer of the 
objective function 



gx_{u^u e ,p e 



m/2 n/2 

E E 

£=-m/2 j=-n/2 



(3.20) 



which quantifies the curvature of the trace. Here Dj? approximates the second derivative in s, and 
u = (u, 0) is defined by 



u m n 



u 



u-t(0) =u ± + u e m o -t(0). 



(3.21) 



Because u± has a weaker effect on the traces than it is more difficult to estimate it. Even if 
there is no noise, the estimation is sensitive to remnants of the traces from the stationary targets 
that have not been perfectly annihilated, the traces from the other moving targets, and the discrete 
sampling of the data. If we took the maximum over I of the sum 



n/2 

E 

j=-n/2 



D r 



as we did in the objective function (3.19), we would seek a saddle point of T to estimate u±. That 
is to say, we would minimize with respect to u± and maximize with respect to £. Our numerical 



simulations show that such an estimation is not robust. To stabilize it, we sum over £ in (3.20) and 
reduce the problem to that of minimizing in u±. 



As an illustration, we present the results for Scene 1. Figure 12 displays Y+ p D r (s,t) using the 



true locations p and velocities u for the two moving targets. Note how after the transformation 
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Figure 12: Left: The traces from the two moving targets in Scene 1. Middle and right: The 
transformed traces T+ p D r (s, t) for the first and second target. 





Figure 13: The objective function g(u) using the traces from the entire Scene 1 (left) and the traces 
from the moving targets in Scene 1 (right). 



the traces become vertical line segments. If we were to search for the range speed u e without 
annihilating the traces of the stationary targets, we would get the objective function g(u) plotted 
on the left in Figure 13 The prominent peak of g is at speed u — 0, because the stationary targets 



dominate the scene. The smaller peaks are due to the two moving targets. When we compute the 
objective function g(u) with the traces from the moving targets alone, we get two large peaks, at 
the range speeds of the targets. The speed in (3.17) is u = ^ e (m o ,0), with u e estimated from the 
first and second peak of g(u), respectively. 

Recall that the cross-range speed u± has a weaker effect on the traces than the range speed, 



it determines only their curvature. We illustrate this in the top right plot of Figure 14, where we 



zoom around the trace of one target, rotated with transformation (3.17) at the estimated range 



speed u e (m o ,0). The bottom row in Figure 14 shows the traces transformed with (3.17), after the 
cross-range speed correction. The trace is now straightened, as seen in the bottom right plot. The 



objective function (3.20) is shown in the middle and right plots of Figure 15, We plot it for the 



two range speeds, the peaks of g shown on the left. The function g± has minima at the cross-range 
speeds of the two targets. 
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Figure 14: Top row: The traces after the travel time transformation (3.17) at the estimated range 



speed of one target. Bottom row: The traces after the travel time transformation (3.17) with the 
estimated cross-range correction of the speed. The right column is a zoom of the left column. 




-40 -20 



Figure 15: Left: The objective function g. Middle and right: The objective function g± for the 
estimated range speeds of the two targets in Scene 1. The minima of g± are the estimates of u±. 



4 Analysis 

In this section we give a brief analysis of our data separation approaches. The analysis is for simple 
scenes, but the methods work in complex scenes, as shown already in the previous section. More 
numerical results are in section [5j We begin with the data model and the scaling regime. Then, we 
analyze the annihilation filters defined in section [3T| The analysis extends to the filters described 



in section [3T4| once we explain how the target velocity affects the traces. We end the section with 
a brief review of the analysis in [2j of the data separation with robust PCA. 
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4.1 Scaling Regime 

We illustrate our scaling regime using the GOTCHA Volumetric SAR data set. The relevant scales 
are the central frequency v Ql the bandwidth £>, the typical range L from the SAR platform to the 
targets, the aperture a, the magnitude |u| of the velocity of the targets, the speed V of the SAR 
platform and the diameter of the imaging set J X . We let L — |r(0) — p G |, and assume that 

B<^v Q , a<L, i? X <L. (4.1) 

We also suppose that the target speed is smaller than that of the SAR platform 

|u| < V. (4.2) 

The GOTCHA parameters satisfy these assumptions: 

• The central frequency is i/q = 9.6GHz and the bandwidth is B — 622MHz. 

• The SAR platform trajectory is circular, at height H = 7.3 km, with radius R = 7.1km 

and thus L w 10km. We consider sub-apertures of one circular degree, which corresponds to 

x 

a = 124m, and image in domains of radius R < 50m. 

• The platform speed is V = 70m/s and the targets move with velocities |u| < 28m/s. 

• The pulse repetition rate is 117 per degree, which means that a pulse is sent every 1.05m, 
and As = 0.015s. 

For a stationary target, we obtain from basic resolution theory [HI [18] that the range resolution 
is c/B = 48cm, and the cross range resolution is XoL/a = 2.5m, with one degree aperture a and 
central wavelength Ao = 3cm. 

4.2 Data model 

We use a simple data model that approximates the scatterers by point targets in the horizontal 
plane and neglects multiple scattering between them. We also suppose that the target trajectories 
can be approximated locally, for the slow time interval defining a single sub-aperture, by the straight 
line segment 

p(s) p(0) + su, 

with constant speed u = (u, 0), and uGl 2 . The traces are modeled by 

N , v 

t} * E { ^\r(T-t(sW fpit ~ (T(S ' MS)) ~ T(S ' ^ ))} ' (4 ' 3) 
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where N is the number of targets. They are located at p q (s) and with reflectivity a q (cj ), for 
q = 1, . . . , N. The model neglects the displacement of the targets during the round trip travel 
time. This is justified in radar applications because the waves travel at the speed of light which is 
several orders of magnitude larger than the speed of the targets. We refer to [3J for a derivation 



of model (4.3) in the scaling regime described above. We simplify it further by approximating the 



amplitude factors 

1 1 



47r|r(s) - p q (s)\ 47rL' 
and assuming that the targets have the same reflectivity 

a q (cj ) = ct(cj ), q = 1,. . . ,N. 

The model of the traces becomes 



D r (s,t)^^^M(s,t), (4.4) 



with 

N 

M{s, t) = Y, fp(t ~ M*, P q (s))), (4.5) 

q=l 



and At given in ( |3.2[ ). Thus, the traces are approximately, up to a multiplicative constant, equal 
to the matrix {M(sj, fy)}, for j = — n/2, . . . , N/2 and £ = —rn/2. . . . ,m/2. We neglect the multi- 
plicative constant hereafter, and call M the matrix of traces. 



We could take any pulse f p in (4.5), but to simplify the calculations in the remainder of the 



section, we assume it Gaussian, modulated by a cosine at the central frequency, 

f p (t) = cos(uj t)e- B2t2 / 2 . (4.6) 

4.3 The annihilation filter 



Recall that the filter Q p defined by (3.4) is a linear operator. Because we model the traces (4.5) 
as a superposition of pulses f p delayed by the travel times to the targets, we can study the effect 
of Qp e on one trace at a time. 

Consider the trace from a stationary or moving target at location Pj(s), for some integer j 
satisfying 1 < j < TV, and denote it by 

M j (s,t) = f p (t-Ar(s,pAs))). (4.7) 
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The transformation defined in (3.1) maps the trace to 

Tf M j ] (8, t) =f p (t + At(s, p e ) - Ar(s, p jS ))) , 
and then we filter it with the difference operator D s , 



(s, t) « A(s, pj(8), p e )f p (t + Ar( S , p e ) - Ar(s, pj( S ))) 



(4.9) 



Here we approximate the finite difference by the derivative in s, denote by f' p the derivative of the 
pulse, and let A be the "annihilation factor" 



d 



A(s,Pj(s),p e ) = J- g [T(8,Pj{s))-T{8,pP 



(4.10) 



Using the scaling assumptions ( |4.1|4.2[ ), we obtain after a straightforward calculation that 
d 



ds 



[t(s,Pj(s)) -T( S ,p e )] 



-u • m e + 



(Vt - u) • FeApj (2Vt - u) • P e u 



+ 



'aVR 1 



(4.11) 



where 



f (0) — p e 

Apj = pj(0) - p e , m e = ,^ n ( ; e| , P e = I - m e mj. 



|r(0) - pt 

The annihilation factor of a stationary target is given by 



A(8,p J ,p*) = -— [ -^ + 0^) 



(4.12) 



It is approximately linear in the cross-range component of the error Ap e of the estimated location 
of the target. For a moving target the factor is 



A(s,pAs),p e 



2V 

c 



u 
V 



Vs U± 

'17 V 



+ 



o( 



au± 

\~cl 



+ o 



aVR 1 
cL 2 



(4.13) 



It is much larger than (4.12) for stationary targets that are in the vicinity of p e , and satisfy the 
condition 



\ApA < L 



(4.14) 



Consequently, the filter Q p defined by (3.4) emphasizes the traces from the moving targets and 



suppresses the traces from the stationary targets near p e , as stated in section 3.1 



Remark 1 In the GOTCHA regime, and for one degree apertures as described in section the 
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cross-range resolution of a preliminary image is approximately 2.5m. Thus, we can assume that 
for the target that we wish to annihilate |Ap^| w 2.5m, and condition (4-H) is satisfied for targets 
moving at range speed that exceeds 0.2ra/s. 



Note also from (4.11), with p e replaced by p G , that the slope of the trace is approximately given 



by the range speed u = u • m , as assumed in section 3.4, as long as 



u > V ~ V— 



If the target is a slower mover, than its trace looks similar to that of a stationary target and it 
cannot be separated. The curvature of the trace is determined by the cross-range velocity of the 



target, as seen from the third term in the right hand side of equation (4.11). 



4.4 Analysis of data separation with robust PCA 

In this section we review briefly from [2] the analysis of the rank of the matrix of traces for simple 
scenes. The goal is to show that indeed, traces from stationary targets should be expected to give a 
low rank contribution whereas those from the moving targets should give a high rank contribution. 



This is the key assumption in the data separation approach described in section 3.2 



4.4.1 One target 

We begin with the case of one target at location p(s) = (p(s),0). It moves at constant speed 
u = (u, 0) over the slow time interval |s| < 5(a), so we can write that 



p(s) = p + su, \s\ < S(a), 

where 

p = p(0). 

The matrix denoted by M is defined by the discrete samples of 

B 2 



M(s, t) — cos [uj (t — Ar(s, p(s)))] exp 



(t-Ar(s,p(s))Y 



(4.15) 



(4.16) 



where Ar is defined by (3.2), and we used the pulse (4.6). The rank of M. is the same as that of 



the symmetric, square matrix y € ]R( n+1 ) x ( n + 1 ), given by 



y = MM\ 



(4.17) 
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Its entries are the discrete samples of the function 



"V 2 i ,00 

y(s,s')= J2 M(s,t q )M(s',t q )*t— dtM(s,t)M(s',t). (4.18) 

/o " — OO 

g=— m/ 2 

We assume here that At is small enough to approximate the Riemann sum over q by the integral 
over t. Since the traces vanish for \t\ > As/2, we can extend the integral to the whole real line. 
We showed in [2] that under our scaling assumptions, the function ( |4.18| ) takes the form 

s/) * c ° s - ^ exp [- ( ^ )2( r* /)2 ] . ( 4 - iq ) 

with dimensionless parameter 

q = 2u - m Q 2Vt -¥ (p- p Q ) | 2u • F Q (p - p Q ) ^ 
c cL cL 



That is to say, the matrix (4.17) is approximately Toeplitz, with diagonals given by the entries in 
the sequence {yj}j e z, 

where 

_(BAs\a\jf_ 

Vj = 2BAt e 4 cos KAs|a|j). (4.21) 

There are many slow time samples in an aperture (n ^> 1), so the matrix 3^ is also large, and we 
can approximate its rank using the asymptotic Szego theory jT9j [6] . 

Note that multiplication of a large Toeplitz matrix with a vector can be interpreted approxi- 
mately as a convolution, which is why the spectrum can be estimated from the symbol 5/(0), the 



Fourier transform of (4.21), 



y(0)= E (4-22) 

j=-oo 

Szego's first limit theorem [6] gives that 

r M(n-fa,fa) 1 f\ t^taWM (aoi\ 

n-S-oo n + 1 Z7T 7_ w lK1 ' K ' iJ 

where l[/3 1; /3 2 ] is the indicator function of the interval [ySi , ^2] and N{n; fa, fa) is the number of 
eigenvalues of ^ that lie in this interval. Using this result, we estimated the rank in [2] as 

rank[^] :=JV(n;e||^|oo,oo), ||y||oo = sup \y{9)\. (4.24) 

0e(-7r,7r) 
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Estimated Rank 
Computed Rank 
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Figure 16: Comparison of the computed and estimated rank of y for a single target. Left: Station- 
ary target at various cross-range positions. Right: Moving target with various velocities. 



Here < e C 1 is a small threshold parameter, and \\y\\oo is of the order of the largest singular 
value of 3^- The Szego theory [THE] gives that this singular value is approximated by the maximum 



of the symbol. The right hand side in (4.23) is calculated explicitly in [2J, and the rank estimate is 

rank [y] 



n + 



r 1 * h f^m^oo) (m)de 

2\a\BAsJ\ogl/e \ 

-LJ v 6 1 , 1 . (4.25) 



mm 



It grows linearly with a, meaning that the larger the offset p — p Q in the cross-range direction is, 
and the faster the target moves, the higher the rank. 



We show in Figure 16 an illustration of this result. The left plot is for a stationary target and 



the right plot is for a moving target. The green line shows the numerically computed rank of the 



matrix y, while the blue line shows its asymptotic estimate (4.25). Note that both curves show 
the same growth. In the case of the stationary target there is a mismatch of the magnitude of the 
computed and estimated rank. This is because n is not large enough in the simulation. A good 
approximation is obtained when yj « for j « n, because then we can approximate the series 



(4.22) that defines the symbol by the truncated sum for indices \j\ < n. This is not the case in this 
simulation, so there is a discrepancy in the estimated rank. However, if we increase the aperture, 
thus increasing n, the result improves. Figure [17] illustrates this fact by showing how the rank 
normalized by the size of the matrix has the predicted asymptote as n increases. 



In the right plot of Figure 16 we compare the computed rank (in green) and the asymptotic 



estimate ( 4.25| ) (in blue) for a moving target located at p = p Q at s = 0. The rank is plotted as a 



function of the range component of the velocity. Note that, as expected, the rank increases with 
the velocity, and therefore with \a\. For the moving target, the asymptotic estimate is in closer 
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Figure 17: Convergence of the rank of y normalized by the size {n + 1). The blue line is the 
computed value and the green line is the asymptotic estimate. The slow time sampling is fixed, so 
the longer the aperture is, the larger n. 

agreement to the computed one. This is because in this case the entries in the sequence {yj}j e z 
decay faster with j. Finally, we can observe that the rank is much larger for the moving target 
than the stationary target, even for small velocities. 



4.4.2 Two targets 



When there are two targets in the scene, we obtain from model (4.5) and the choice (4.6) of the 
pulse that 



M{s,t) = ^cos [uo (t — Ar(s,Pj))] exp 



-^(t-Ar^p,)) 2 



(4.26) 



We restrict the discussion to the case of two stationary targets. The extension to moving targets 



is straightforward. It amounts to adjusting the parameters a± and a 2 defined below in (4.28), by 



adding two linear terms in the target velocity, as in equation (4.20). 



We obtain after a calculation given in detail in [2J that the entries of the matrix (4.17) are the 
discrete slow time samples of the function 



y(s,s') 



7T 



2BAt 



cos[uj aj(s — s')] exp 



3=1 



(Baj) 2 (s - s' 



A21 



+ cos[oj (ais — a2S + /?)] exp 
+ cos[uj {ais — a2S + /?)] exp 



B 2 (a 1 s-a 2 s f + /3) 2 



B 2 (a 1 s f - a 2 s + 0) 
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(4.27) 
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where 



and 



2Vt • KiPj - p Q ) . 
= ^ , J = 1,2, 



o 2 A f ^ [m G • (0. - pj] 2 



2L 



(4.28) 



(4.29) 



Thus, the matrix y has a more complicated structure. The first term in (4.27), the sum over j 



gives rise to a Toeplitz matrix, as before. The last two terms give matrices that are approximately 
g-Toeplitz or g-Hankel, depending on the sign of the ratio g^/cki. 

Definition 1 An (n + 1) x (n + 1) g-Hankel matrix % with shift g E Z + is defined by a sequence 
{hj}jeN ^ 

Ujt = j,£ = 1, . . . ,n + 1. 

T/ie matrix is Hankel when g — 1. A g-Toeplitz matrix is defined similarly, by replacing g with —g. 



It is easier to analyze the spectrum of the matrix y, when the reference point p Q is chosen so 



that 



Then, y is given by 



OL2 



< 0, and g :- 



OL2 



E N. 



y = T + H + U T , 



with the Toeplitz and g-Hankel matrices 

7j£ = yj-£, Hjt = j,£ = 1, . . . ,n+ 1. 

The sequences {yj}j^z and {^}jGN that define these matrices are 

2 

22JAt ^ C ° S ( £J ° As l a ^lj) e 



(4.30) 
(4.31) 

(4.32) 



(BA S |a^|j)^ 



£=1 



and 



7T [BAs\ ai \(j+C)] z 

e 4 cos[cj A5|ai|(j + QJ, C 



If the range offsets from the reference point p Q are sufficiently large, so that 

2^(-l) l m o • \Pt ~ Po) + 2L 



(4.33) 



(4.34) 
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»1, 



then the g-Hankel matrix has small entries, and y is approximately Toeplitz, as in the single target 
case. The meaning of this is that the difference of the travel times between the SAR platform and 
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Figure 18: Computed and estimated rank of matrix y for a scene with two stationary targets. On 
the left, one target is fixed at location p ± = (5, 5, 0)m and the location of the other varies on the 
line segment between (— 5,0.01,0)m and (—5,30, 0)m. On the right, one target is fixed at location 
Py — (0.15, 5, 0)m and the location of the other varies on the line segment between (—0.15, 0.01, 0)m 
and (-0.15,30,0)m. 



such targets is larger than the compressed pulse width, which makes their interaction in (4.27) 
negligible. However, if the range offsets are small, the interaction plays a role and the structure of 
the matrix y is given by equation (4.31). In this case, the estimate of the rank of y follows from 
the recent results in [211 UHl El- They state that the g-Hankel terms T~L + H T have a negligible 
effect on the rank in the limit n — » oo. See [2j for more details. Thus, we can conclude in either 
case that the rank estimate of y is given by equation (4.23), in terms of the symbol y{9) defined 
by (4.22), using the sequence {yj}j^z with entries (4.33). 

We plot in Figure 18 the computed and estimated rank for two different configurations of the 
stationary targets. The results on the left consider one target at location p x — (5, 5, 0)m. The 
second target location is varied between (— 5,0.01,0)m and (—5,30, 0)m. The range separation is 
large, equal to 10m, so the entries of the g-Hankel matrix H are negligible because B\f3\ = 41.47 ^> 1. 
The plot on the right considers one target at p r = (0.15, 5, 0)m and a second target at location the 
p 2 varying between (—0.15, 0.01, 0)m and (—0.15, 30, 0) m. The range separation is now small, and 
the g-Hankel matrix % is no longer negligible. Nevertheless, as predicted by the asymptotic theory, 
the rank of matrix y behaves essentially the same as before. 



4.4.3 Conclusions of the analysis of rank 

Several conclusions can be made from the analysis summarized above. First, the traces from a 
moving target form a high rank matrix. For a fixed number n of slow time samples, the rank grows 
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linearly with the target speed, until the matrix becomes full rank. The traces from a stationary 
target form a low rank matrix if the scene is small enough. The rank grows linearly with the 
cross-range offset of the target from the reference point p G , but the growth is at a slower rate than 
for a moving target. 

We also analyzed the rank for scenes with two stationary targets and found that it increases 
linearly with the cross-range offsets of the targets from p Q . The expectation, confirmed by the 
numerical simulations, is that the rank increases with the number of targets. This is why the data 
separation with robust PCA should be done in successive small time windows, with each window 
containing the traces from only a few stationary targets. These traces give a matrix that is low 
rank, and thus can be separated from the traces due to moving targets. The results shown in Figure 
[9] illustrate this point. 



5 Numerical Simulations 



We present here more numerical simulations that illustrate the performance of the annihilation 



filters defined in sections |3.1| and |3.4| We have presented already in sections |3.2| and |3.3| three 
different simulations for data separation with robust PCA. More results can be found in [2]. 



5.1 Setup of the simulations 



We use the GOTCHA Volumetric SAR setup as described in section 4.1 The data traces are 



generated with the model (4.3). In all the simulations, the point targets are assumed to have 
identical reflectivity cr q = 1. The motion estimation results are obtained with the phase space 
algorithm introduced in [3J. This algorithm requires that we know the location of the target at one 
moment during the measurement window. We choose it at the center of the aperture, which is why 
there is no error in the target trajectory at s = 0. 

For the robust PCA results, the principal component pursuit optimization is solved with an 
augmented Lagrangian approach. It requires the computation of the top few singular values and 
corresponding singular vectors of large and sparse matrices, which we do with the software package 
PROPACK. 



5.2 Annihilation Filters 



We give two examples of data separation with the annihilation filters described in section 3.1 In 
the first example, we have a complex scene composed of ten stationary scatterers in an imaging 
region of 50x50m 2 , and a moving target with velocity u = 28/\/2(l 5 1, 0)m/s. We plot in the top 
row of Figure 19 the data traces, and in the bottom row the filtered traces. The right column is the 
same as the left except it is plotted in decibel (dB) scale to increase the contrast of the figure. The 
second example considers a scene with thirty stationary targets and the same one moving target. 
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Figure 19: Data separation with the annihilation filter for a scene with ten stationary targets and 
a moving one. The top row shows the data traces and the bottom row the trace from the moving 
target, as separated by the algorithm. The right column is the left column plotted in dB scale. 



The results are in Figure 20, In both examples we used the exact stationary scatterer locations in 
the definition of the filters (|3.4[). The errors in the estimated target trajectories, with the separated 



trace from the moving target for each example, are shown in Figure 21 



Note that the filtering and the motion estimation results are worse in the second example, 
where there are more stationary targets to remove, by applying repeatedly filters like (3.4) to the 
traces. Computing many times numerical approximations of slow time derivatives causes errors to 
accumulate. This is specially because we have high frequency signals, and the sampling in slow 
time is limited by the pulse repetition frequency of the SAR system. 



5.3 Data separation for Scene 1 



We already saw in section 3.3, Figure 11, the separation of the traces from the twenty stationary 
and two moving targets in Scene 1. Here we show the further separation of the traces from the 
moving targets, as described in section |3~l} The first moving target is located at (0, 0, 0)m at s = 0, 



and is moving in the horizontal plane with velocity u = ^(1, l)m/s. The second target is located 
at (—5, 5, 0)m at s = 0, and moves with velocity u = ^| (—1, \/2)m/s. 

We apply the method described in section |3.3| to the traces shown in the right plot of Figure 
[TT| separated by robust PCA. We estimate the velocity of the first target, apply the travel time 



transformation (3.17) to the traces, and then use the robust PCA approach to do the separation. 
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Figure 20: Data separation with the annihilation filter for a scene with thirty stationary targets 
and a moving one. The top row shows the data traces and the bottom row the trace from the 
moving target, as separated by the algorithm. The right column is the left column plotted in dB 
scale. 




Figure 21: Error of the estimated target trajectory in a complex scene with ten stationary targets 
(left) and thirty stationary targets (right). The motion estimation is done with the filtered data 
shown in the bottom row of Figures [19] and [20j 



The output of the data separation is shown in Figure 22 The traces from the stationary targets, 
as obtained by the robust PC A method are on the left. The traces from the two moving targets 



are in the middle and right plots. The images obtained with these traces are in Figure 23 All the 
targets, stationary and moving are now in focus. The peak of the image for the second moving 
target is slightly displaced in cross-range due to the error in its velocity estimation after separation. 
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Fast time (seconds) xio -7 Fast time (seconds) xio -7 Fast time (seconds) xio -7 

Figure 22: Data separation for Scene 1. Left: traces from the stationary targets. Middle and right: 
traces from the moving targets. 




Range (meters) Range (meters) Range (meters) 



Figure 23: Images for Scene 1. Left: image with the traces from the stationary targets. Middle 
and right: image with the traces the moving targets. 

6 Summary 

We presented two approaches for separating the pulse and range compressed echoes from stationary 
and moving targets in a complex scene using a single synthetic aperture radar (SAR) antenna. The 
separation is required for estimating the motion of the targets and focusing the images. 

The first approach subtracts the echoes from the stationary targets one by one, by performing 
the following four steps: (1) It estimates the stationary target location p e from a preliminary image. 

-*e 

(2) It transforms the data using a map T:j_ so that the echo from the target becomes approximately 
independent of the SAR antenna location on the flight path. (3) It annihilates the target echo by 
exploiting this independence. (4) It undoes the transformation at step (2) by applying the map 
to the remaining echoes. 

The second approach is based on the robust principle component analysis (PC A) method, which 
is designed to decompose a matrix Ai into its low rank part and high rank but sparse part. In our 
context, the entries of the matrix A4 are the samples of the pulse and range compressed echoes 
measured by the SAR antenna. The main observation is that with appropriate pre-processing and 
windowing, the discrete samples of the stationary target echoes form a low rank matrix, while those 
from a few moving targets form a high rank sparse matrix. Thus, they can be separated with robust 
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PCA. 

We presented a brief analysis of the two methods, and a numerical comparison of their perfor- 
mance. Each method has its advantages and disadvantages, but they complement each other. Thus, 
we can combine them to improve the data separation results for complex and extended imaging 
scenes. 
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